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Abstract:  Cazenovia  Creek,  in  Western  New  York,  is 
the  largest  tributary  of  the  Buffalo  River.  Breakup  ice 
jams  form  along  the  lower  basin  nearly  every  year  dur¬ 
ing  mid-winter  or  spring  thaws,  and  ice-jam  flooding 
occurs  in  the  City  of  Buffalo  and  the  Town  of  West  Sen¬ 
eca  about  every  2-3  years.  This  report  describes  phys¬ 
ical  model  tests  and  design  recommendations  for  a  new 
ice-control  structure  (ICS)  for  Cazenovia  Creek.  The 
recommended  structure  consists  of  nine  10-ft-tall  x  5-ft- 
diameter  cylindrical  piers  spaced  across  the  main  chan¬ 
nel,  and  it  uses  the  adjoining  treed  floodplain  as  a  natu¬ 


ral  bypass  channel.  Also  described  are  results  from  a 
numerical  ice-hydraulic  model  to  determine  the  extent 
of  flooding  induced  upstream  of  the  new  ICS.  Although 
few  structures  are  affected,  the  ice  jam  held  by  the  ICS 
will  cause  minor  flooding  of  properties  abutting  the  creek. 
However,  the  stream-wise  extent  of  this  flooding  will 
decrease  during  an  event  as  melting  and  washouts 
reduce  the  volume  of  ice  in  the  jam.  The  structure  bal¬ 
ances  the  need  to  protect  downstream  areas  from  nat¬ 
ural  ice-jam  flooding  and  the  need  to  minimize  upstream 
flooding  induced  by  the  retained  ice. 
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CONVERSION  FACTORS:  U.S.  CUSTOWARY  TO  METRIC 
(SI)  UNITS  OF  MEASUREMENT 

These  conversion  factors  include  all  the  significant  digits  given  in  the 
conversion  tables  in  the  ASTM  Metric  Practice  Guide  (E  380),  which 
has  been  approved  for  use  by  the  Department  of  Defense.  Converted 
values  should  be  roimded  to  have  the  same  precision  as  the  original 
(see  E  380). 


Multiply 

By 
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inch 
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foots 
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mile 
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joule/kilogram  kelvin 

Btu/ftS  °F 
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degrees  Fahrenheit 

fc  =  0F-”32)/1.8 

degrees  Celsius 
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INTRODUCTION 


Project  Overview 

Cazenovia  Creek,  located  in 
western  New  York  State,  is  the  larg¬ 
est  tributary  of  the  Buffalo  River 
(Fig.  1 ).  Its  basin  is  subject  to  heavy 
lake-effect  snow  from  Lake  Erie, 
and  spring  runoff  typically  has  a 
strong  snowmelt  component.  Dur¬ 
ing  most  winters,  ice  covers  80- 
100%  of  the  surface  of  the  creek, 
and  its  average  thickness  ranges 
from  1  to  2  ft.  Breakup  ice  jams 
form  nearly  every  year  during  mid¬ 
winter  or  spring  thaws,  and  flood¬ 
ing  occurs  along  the  lower  basin  in 
the  City  of  Buffalo  and  the  Town  of 
West  Seneca  about  every  2-3  years. 
Some  of  these  floods  have  caused 
damages  exceeding  $1,000,000,  and 
one  death  has  occurred. 

During  the  1980’s,  the  U.S. 
Army  Engineer  District,  Buffalo 
(CELRB),  investigated  the  feasibil¬ 
ity  of  building  an  ice-control  struc¬ 
ture  (ICS)  to  reduce  ice-jam  flood 
damages  along  Cazenovia  Creek  in 
the  Town  of  West  Seneca  (U.S. 
Army  1986a,  1986b).  Based  on 
model  tests  at  CRREL  (Gooch  and 
Deck  1990),  the  ICS  chosen  con¬ 
sisted  of  a  6-ft-high  x  250-ft-wide 


Figure  1.  Cazenovia  Creek  drainage  basin.  (From  USACE  1975.)  The  new 
cylindrical-pier  ICS  would  be  located  at  the  same  site  originally  proposed  for  a 
welr-with-plers  ICS  or  “Ice  retention  structure”  (USACE  1986a,  1986b). 
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Figure  2.  Concept  drawing  of  original  Cazenovia  Creek  !CS.  Dimensions  and  elevations  shown 
are  prototype  feet.  Not  shown  is  a  low-flow  gated  opening  in  the  weir  to  allow  fish  passage  from  spring  to 
early  winter.  In  its  final  design,  the  ICS  piers  were  12  ft  long.  A  pool  for  storage  of  ice  floes  was  to  be 
excavated  pool  at  elevation  635;  it  was  400  ft  wide  at  the  ICS  and  tapered  linearly  to  the  natural  channel 
width  (approximately  160  ft)  at  a  distance  of  600  ft  upstream.  The  project  also  required  extensive  grading 
and  riprap  protection  for  the  floodway  bypass  channel. 


weir  with  nine  ice-retaining  piers,  an  excavated  pool  to 
store  ice  pieces,  and  a  prepared  floodway  to  bypass  the 
peak  flow  (Fig,  2).  The  project  was  estimated  to  cost 
$2,100,000  (1987  dollars),  reduce  average  annual  ice- 
jam  flood  damages  by  $320,000  (1987  dollars),  and 
achieve  a  benefit/cost  ratio  of  about  1 .5  over  the  life  of 
the  structure.  Nevertheless,  it  was  not  built.  Federal  law 
was  changed  at  that  time  to  require  25%  cost  sharing 
by  the  local  sponsor.  The  Town  of  West  Seneca  could 
not  afford  this  amount  on  its  own,  and  New  York  State 
had  not  yet  modified  its  laws  to  participate  as  a  spon¬ 
sor. 

Many  communities  along  small  northern  rivers 
experience  severe  ice-jam  flooding.  Although  flood  dam¬ 
ages  can  be  significant,  ICS’s  that  include  weirs,  dams, 
or  levees  are  often  too  expensive  to  achieve  benefit/cost 
ratios  over  unity.  In  1991,  CREEL  began  research  spe¬ 
cifically  to  reduce  ICS  costs  for  small  rivers.  The  result 
for  breakup  ice  jams  was  the  sloped-block  ICS,  devel¬ 
oped  using  model  tests,  and  demonstrated  in  the  Lam¬ 
oille  River  in  Hardwick,  Vermont  (Lever  et  al.  1997, 
Lever  and  Gooch  1998).  The  Hardwick  ICS  consists  of 
four  massive,  sloped  blocks  placed  across  the  river  ad¬ 
jacent  to  a  treed  floodplain.  The  blocks  arrest  a  break¬ 
up  ice  run  and  form  a  stable,  partially  grounded  ice  jam. 
Trees  on  the  floodplain  retain  ice  pieces  in  the  river 
channel  while  allowing  flow  to  bypass  the  structure. 
The  large  gaps  between  blocks  (14  ft)  allow  easy  fish 
and  canoe  passage. 

The  sloped-block  ICS  has  performed  well,  and  no 
ice-jam  flooding  has  occurred  in  Hardwick  since  its 
construction  in  1994.  However,  for  ice  less  than  1  ft 
thick,  the  ICS  released  the  ice  jams  after  several  hours. 
This  does  not  pose  a  problem  in  Hardwick  and  is  a  con¬ 
sequence  of  the  safety-valve  feature  of  the  sloped  blocks 
to  minimize  peak  loads  (and  construction  costs).  But 
many  areas  downstream  of  the  Cazenovia  ICS  site,  in 


both  West  Seneca  and  Buffalo,  are  heavily  developed. 
Consequently,  we  sought  greater  ice-retaining  capac¬ 
ity  for  the  Cazenovia  ICS. 

Figure  3  shows  the  new  ICS  concept  proposed  for 
Cazenovia  Creek.  It  consists  of  evenly  spaced  cylin¬ 
drical  piers  anchored  in  main  channel  and  does  not 
include  a  weir.  The  adjacent  treed  floodplain  is  left  intact 
to  act  as  a  flow-bypass  channel.  In  proposing  this  ICS, 
we  expected  that  vertical  piers  spaced  the  same  dis¬ 
tance  apart  as  the  Hardwick  ICS  blocks  would  offer 
significantly  better  ice-retention  capability,  especially 
for  thinner  ice.  Nevertheless,  we  needed  to  quantify  this 
and  verify  that  piers  alone  would  offer  the  same  or  bet¬ 
ter  performance  as  the  original  weir-with-piers  ICS. 
Also,  we  needed  to  determine  design  ice  loads  on  the 
piers.  We  therefore  conducted  physical  model  tests  of 
the  new  ICS. 


Figure  3.  Concept  drawing  of  new  Cazenovia  Creek 
ICS.  Cylindrical  piers  of  5  ft  diameter  X 1 0  ft  tall  with  1 0-  to 
14-ft  gaps  anchored  in  main  channel  (no  weir)  arrest  a 
breakup  ice  run  and  retain  the  ice  jam  throughout  an  event. 
Treed  floodplain  bypasses  flood  flow  but  not  ice  floes. 
Raised  riprap  embankment  near  ICS  prevents  scour  and 
ensures  floodplain  flow  returns  to  main  channel  down¬ 
stream  of  ICS.  The  riverbed  at  the  ICS  site  is  bedrock. 
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These  test  results  showed  that,  as  expected,  the  cylin¬ 
drical-pier  ICS  will  retain  an  ice  jam  at  discharges  much 
higher  than  that  needed  to  release  a  natural  ice  jam. 
Consequently,  upstream  water  levels  will  increase  com¬ 
pared  to  existing  conditions.  We  used  a  numerical  ico- 
hydraulic  model  to  determine  the  maximum  water  sur¬ 
face  elevations  upstream  of  the  structure  and  the  loca¬ 
tion  beyond  which  the  ICS  should  cause  no  effect. 
CELRB  then  used  these  results  to  establish  real  estate 
requirements  for  project  construction. 

This  report  describes  the  model  tests  of  the  new 
Cazenovia  Creek  ICS,  the  design  recommendations 
based  on  those  tests,  and  the  ICS’s  expected  upstream 
effects.  Where  possible,  we  compare  results  with  the 
original  weir-with-piers  ICS.  Note  that  all  physical  units 
reported  here  are  prototype  scale,  including  those 
derived  from  model  tests.  All  elevations  quoted  are 
National  Geodetic  Vertical  Datum  (NGVD). 

Ice-hydraulic  conditions 

Figure  1  shows  a  map  of  the  Cazenovia  Creek  water¬ 
shed  in  western  New  York.  The  drainage  area  is  144 
mi^.  The  main  stem  of  the  Cazenovia  Creek  is  about  1 7 
miles  long  from  the  confluence  of  the  creek’s  east  and 
west  branches  in  East  Aurora  to  the  Buffalo  River.  Its 
average  slope  along  this  section  is  0.0026  or  about  14 
ft/mile. 

As  mentioned  earlier,  snowmelt  and  rainfall  events 
in  mid-winter  and  early  spring  have  caused  numerous 
ice-jam  floods  along  Cazenovia  Creek  in  West  Seneca 
and  Buffalo.  The  main  stem  contributes  most  of  the  ice 
volume  to  these  jams.  Increasing  discharge  fractures 
the  ice  cover,  lifts  it  free  of  the  banks,  and  sends  it  down¬ 
stream  as  a  surge  of  water  and  ice  floes.  Shear  walls 
(i.e.,  walls  of  packed  ice)  have  been  observed  from  Tran¬ 
sit  Road  to  Mill  Road,  indicating  at  least  temporary 
formation  of  ice  jams.  However,  the  damage-causing 
ice  jams  form  when  the  ice  run  encounters  the  strong, 
thick  ice  sheet  along  the  flatter  reach  below  Mill  Road. 
These  jams  may  release  after  several  hours  and  reform 
downstream  until  discharge  is  sufficient  to  clear  ice  into 
the  Buffalo  River  and  thence  to  Lake  Erie. 

About  134  mi^  of  drainage  area  lies  upstream  of  the 
USGS  Ebenezer  gage  at  Ridge  Road.  Figure  4  shows 
stage  hydrographs  from  the  gage  for  ice  jam  floods 
occurring  on  2  March  1972  and  23  February  1985.  The 
hydrographs  show  characteristics  typical  of  ice-jam 
events.  The  ice-affected  stage  increases  slowly  as  the 
event  begins.  When  an  ice  run  passes  the  gage,  the 
measured  stage  rises  sharply  and  falls  abruptly.  The  gage 
drops  back  to  its  open-water  rating  curve  after  the  ice 
moves  far  downstream.  The  remainder  of  the  hydro¬ 
graph,  including  the  broad  peak,  is  unaffected  by  ice 
and  can  provide  accurate  discharge  measurements.  For 


a.  2  March  1972.  (After  U.S.  Army  1972.) 


February  19B5 

b.  23  February  1985 

Figure  4.  Hydrographs  from  USGS  Ebenezer  gage 
for  two  ice-jam  events. 

the  1972  and  1985  events,  the  discharges  correspond¬ 
ing  to  the  broad,  open-water  peaks  were  about  9500 
and  9000  cfs,  respectively.  The  maximum  discharge  of 
record  (since  1941)  was  13,500  cfs  on  1  March  1955. 
The  maximum  discharge  occurring  during  a  known  ice- 
jam  event  was  12,600  cfs  on  22  January  1959.  Ice  thick¬ 
ness  prior  to  breakup  ranges  widely,  but  few  flooding 
problems  have  occurred  for  ice  thinner  than  about  10 
in. 

The  site  selected  for  both  ICS  concepts  lies  about 
2300  ft  upstream  of  Mill  Road  (Fig.  5).  The  creek  at 
this  location  drains  129  mP,  its  main  channel  is  about 
150  ft  wide  and  7  ft  deep,  and  the  imdeveloped,  treed 
floodplain  is  about  400  ft  wide.  The  channel  bed  through 
the  2000-ft-long  ICS  reach  consists  of  shale  bedrock 
and  has  an  average  slope  of  0.001;  this  slope  increases 
to  about  0.003  for  the  next  4  miles  upstream.  For  1 1 
miles  upstream  of  Mill  Road,  the  creek  flows  through 
an  incised  valley  that  contains  very  little  developed 
property. 
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Figure  5.  Lower  basin  of  Cazenovra  Creek  and  the  site  of  the  ice-control  structures  discussed  here. 


iCS  DESIGN 

Original  ICS  features  and  function 

The  original  weir-with-piers  ICS  (see  Fig.  2)  was 
based  on  a  similar  structure  installed  on  the  Ste.  Anne 
River  in  St.  Raymond,  Quebec  (Deck  1984).  In  early 
winter,  a  low-flow  opening  would  be  closed  to  create  a 
pool  and  form  an  ice  sheet.  The  weir  and  excavated 
pool  were  sized  to  reduce  water  velocities  sufficiently 
to  store  ice  floes  arriving  from  upstream  during  break¬ 
up.  Piers  protruded  above  the  weir  to  stabilize  the  ice 
sheet  and  prevent  passage  of  ice  floes  over  the  weir.  A 
prepared  floodway  helped  reduce  velocities  over  the 
weir  at  high  flows. 

Gooch  and  Deck  (1990)  conducted  physical  model 
tests  to  optimize  this  ICS.  They  determined  that  ice 
breakup  typically  begins  at  about  1500  cfs  and  that  the 
seven  ice-jam  floods  on  Cazenovia  Creek  from  1971- 
82  had  peak  discharges  less  than  6000  cfs.  The  Corps 
accepted  this  value  as  the  design  discharge  for  the  orig¬ 
inal  ICS  (U.S.  Army  1986a),  although  the  hydrographs 
for  the  1972  and  1985  events  indicate  that  discharge 
can  exceed  this  value  during  later  stages  of  breakup 
events  (see  Fig.  4). 

Gooch  and  Deck  modeled  about  4200  ft  of  river  us¬ 
ing  scales  of  1 :40  horizontal  and  1:10  vertical,  and  used 
urea-doped  ice  to  scale  ice  flexural  strength  (prototype 
target  800  kPa  or  120  psi).  Test  results  indicated  that  a 
6-ft-high  weir  with  piers  retained  ice  more  effectively 
than  an  8-ft-high  weir  with  no  piers  (U.S.  Army  1986a). 
For  fewer  than  five  piers,  the  ice  sheet  would  extrude 
past  the  piers,  or  break  up  in  front  of  them,  and  the  ice 
floes  accumulated  in  the  pool  would  spill  over  the  weir 
(Gooch  and  Deck  1990).  Even  with  no  floodway,  the 
6-ft-high  weir  with  five  or  nine  piers  was  found  to  retain 


ice  to  the  design  discharge  (6000  cfs,  which  was  also 
the  maximum  discharge  possible  in  the  model).  Flood¬ 
way  flow  reduced  the  flow  over  the  weir  by  about  20% 
at  the  design  discharge,  providing  some  safety  margin. 
For  its  final  design,  CELRB  selected  the  6-ft-high  weir 
with  nine  piers  and  the  prepared  floodway,  reasoning 
that  increased  reliability  outweighed  the  small  incre¬ 
mental  cost  (U.S.  Army  1986a).  Existing  trees  were  to 
be  left  intact  along  the  floodway  to  prevent  ice  passage 
downstream  at  high  flows. 

The  weir-with-piers  ICS  would  have  required  annual 
maintenance.  All  the  estimated  bed-load  sediment  and 
some  of  the  suspended  sediment  were  expected  to 
deposit  in  the  pool  upstream  of  the  ICS.  Dredging  of 
about  4000  ft^  of  sediment  annually  would  have  been 
needed  to  ensure  adequate  ice-storage  volume  (U.S. 
Army  1986a). 

With  the  inclusion  of  a  low-flow  opening,  no  signif¬ 
icant  long-term  environmental  impacts  on  water  quality 
or  stream  ecosystems  were  expected  to  result  from  the 
original  ICS.  Also,  open-water  flooding  upstream  was 
not  expected  to  increase  because  increased  water- 
surface  elevations  caused  by  the  structure  would  extend 
no  further  than  4500  ft  upstream  (U.S.  Army  1986a). 
Upstream  effects  with  ice  included  were  not  examined, 
based  on  the  assumption  that  all  upstream  ice  would 
collect  within  the  ICS  pool. 

Physical  model  of  cylindrical-pier  ICS 

The  main  requirement  for  the  new,  cylindrical-pier 
ICS  is  to  perform  as  well  as  the  weir-with-pier  ICS. 
That  is,  it  should  arrest  breakup  ice  runs  and  retain  the 
resulting  ice  jams  (for  ice  thicker  than  about  0.8  ft)  for 
the  6000-cfs  design  discharge  (12-hour  rise  time).  In 
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addition,  any  ice  releases  should  be  gradual,  not  cata¬ 
strophic,  and  the  ICS  should  have  minimal  environmen¬ 
tal  impact  and  low  operation  and  maintenance  require¬ 
ments.  The  objectives  of  the  model  tests  were  to  quan¬ 
tify  the  ice-retention  capability  of  the  new  ICS,  opti¬ 
mize  pier  spacing,  and  determine  design  ice  loads. 

We  designed  the  physical  model  based  on  our  expe¬ 
rience  with  the  sloped-block  ICS  (Lever  et  al.  1997). 
In  particular,  we  selected  a  large  scale-factor  (1:15  un¬ 
distorted)  to  study  the  unsteady  ice-jam  processes  that 
occur  near  the  ICS.  For  vertical  piers,  the  gap  between 
piers,  rather  than  their  diameter  or  shape,  should  have 
the  strongest  effect  on  the  ICS’s  ice-retention  capabili¬ 
ty.  Therefore,  we  tested  only  5-ft-diameter  cylindrical 
piers  and  varied  gap  width  from  10-14  ft.  Tests  of  the 
sloped-block  ICS  indicated  that  they  could  retain  ice 
jams  that  rose  about  3  ft  above  the  tops  of  the  blocks. 
We  therefore  selected  a  single  pier  height,  10  ft  above 
the  average  bed  elevation  at  the  ICS  cross-section, 
which  placed  the  top  of  the  piers  about  3  ft  above  the 


adjacent  floodplain  elevation.  We  also  tested  8-ft-talI  x 
5-ft-wide,  45°  sloped  blocks  to  compare  their  perfor¬ 
mance  with  the  cylindrical  piers  and  to  compare  their 
model  behavior  with  the  field  behavior  of  the  Hard¬ 
wick  ICS. 

Figure  6  shows  a  plan  view  of  the  1 : 1 5-scale  refrig¬ 
erated  hydraulic  model  of  Cazenovia  Creek  used  to 
develop  the  cylindrical-pier  ICS.  The  model  covered 
about  2000  ft  of  river  length  and  560  ft  of  width,  and 
included  most  of  the  adjoining  floodplain.  The  ICS  was 
located  in  a  channel  section  about  150  ft  wide  (Fig.  7). 
Time  series  measurements  in  the  model  included  inflow, 
floodplain,  and  channel  discharges,  channel  and  flood- 
plain  water-surface  elevations,  and  loads  on  five  piers. 
Sample  rate  was  2.6  Hz  (prototype  scale).  Single-axis 
load  cells  measured  downstream  overturning  moments 
on  four  of  the  ICS  piers  (M1-M4),  while  a  six-axis  load 
cell  measured  downstream  and  transverse  overturning 
moments  and  downstream  force  on  one  pier  (Fig.  7). 
We  calibrated  the  load  cells  before  each  test  and  zeroed 
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Figure  6.  Plan  view  of  1:1 5-scale  refrigerated  hydraulic  model  of  Cazenovia  Creek.  Chan¬ 
nel  width  at  the  ICS  is  150  ft,  and  the  overali  iength  modeied  is  about  2000  ft. 


Figure  7.  Elevation  view  of  ICS  consisting  of  nine  5-ft-dlameter  x  10-ft-tall 
cylindrical  piers  with  12-ft  gaps.  Dimensions  and  elevations  (EL)  are  prototype 
feet.  Piers  M1-M4  were  equipped  with  single-axis  load  ceils  to  measure  downstream 
overturning  moments,  and  a  six-axis  load  ceil  to  measured  downstream  and  trans¬ 
verse  moments  and  downstream  force. 
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Table  1.  Froude  scaling  expres- 

sions  between  model  (m)  and 

prototype  (p). 

Parameter 

Scaling  expression 

length 

Lp/L„  =  15 

velocity 

Vp/V,n  =  151/2 

time 

Wm  =  151/2 

discharge 

Qp/Qn,=  155/2 

force 

Fp//=m  =  153 

moment 

=  IS-* 

ice  strength 

them  to  remove  the  deadweight  of  the  piers.  Although 
they  were  fewer  than  the  number  of  piers,  we  distributed 
these  load  cells  to  obtain  representative  results.  We  also 
manually  measured  ice  thickness,  flexural  strength,  and 
piece-size  distribution,  and  video  cameras  recorded  the 
tests.  We  used  Froude  scaling  to  relate  model  and  pro¬ 
totype  values  (Table  1). 

Most  tests  were  conducted  to  simulate  ice  breakup 
along  Cazenovia  Creek.  Ice  covers  about  1-2  ft  thick 
were  grown  on  the  river  channel  and  in  the  ice-supply 
flume.  We  would  manually  break  the  latter  ice  cover 
and  feed  it  into  the  head  of  the  model  behind  a  porous 
gate.  We  would  also  partially  break  the  ice  cover  in  the 
river  channel,  usually  leaving  about  160  ft  intact  imme¬ 
diately  upstream  of  the  ICS.  After  establishing  a  base 
flow  of  1000-1500  cfs,  we  would  release  the  gate  to 
simulate  arrival  of  an  ice  run  at  the  ICS  reach.  This  ice 
run  would  form  an  ice  jam,  which  would  then  collapse 
in  stages  onto  the  ICS  as  we  increased  discharge.  The 
tests  continued  until  the  ice  jam  released  or  we  reached 
maximum  discharge  (less  than  8000  cfs)  after  about 
10-12  hours  (prototype).  Figure  8  shows  a  typical  ice 
jam  formed  by  the  cylindrical-pier  ICS  and  the  result¬ 
ing  floodplain  flow. 

Unless  ice  conditions  are  very  mild  (and  thus  pose 


little  ice-jam  threat),  an  ice  sheet  will  be  present  up¬ 
stream  of  the  prototype  ICS  at  breakup.  Nevertheless, 
we  conducted  a  few  tests  with  no  ice  cover  in  the  main 
channel  upstream  of  the  ICS  because  they  represent  a 
simpler  case  for  validating  a  numerical  model  of  break¬ 
up  jamming  at  an  ICS  (Hopkins  et  al.  1996). 

Model  results 

Table  2  summarizes  the  test  conditions  and  results 
of  the  physical  model  study.  We  observed  three  main 
outcomes  for  tests  with  ice  covers  initially  upstream  of 
the  ICS:  1)  no  release  prior  to  the  maximum  pumping 
capacity,  2)  slow  washout  of  ice  floes  through  a  gap,  and 
3)  rapid  release  of  most  of  the  jam  through  more  than  one 
gap.  For  discharges  below  6000  cfs  (the  design  discharge), 
there  were  no  releases  through  the  cylindrical-pier  ICS, 
regardless  of  gap  width.  Releases  took  place  above  the 
design  discharge  on  only  three  occasions  (tests  18,  20, 
and  23),  but  in  the  only  test  where  ice  released  over  the 
cylindrical  piers  (test  18),  the  ice  was  relatively  thin 
(0.89  ft)  and  the  discharge  had  reached  9100  cfs.  Wash¬ 
outs  were  more  common  and  occurred  for  all  three  gaps 
tested,  and  it  is  difficult  to  discern  a  quantitative  differ¬ 
ence  in  performance.  Qualitatively,  however,  the  14-ft- 
gap  cylindrical-pier  ICS  tended  to  produce  brief  wash¬ 
outs  at  lower  discharges  and  longer-duration  washouts 
near  maximum  discharge,  especially  for  thinner  or 
weaker  ice.  We  saw  no  qualitative  or  quantitative  dif¬ 
ference  between  the  performance  of  the  cylindrical-pier 
ICS  for  the  two  smaller  gaps. 

Comparable  tests  of  the  sloped-block  and  cylindrical- 
pier  ICS’s  at  the  same  pier  gap  (14  ft)  illustrate  the  in¬ 
creased  ice-reteqtion  capability  of  the  latter  structure. 
For  moderate  ice  thickness  (1.0-1 . 1  ft)  and  strength  (70- 
100  psi),  the  ice  jams  washed  out  between  or  released 
over  the  sloped-blocks  at  discharges  below  the  design 
value  (tests  13  and  14)  but  were  held  well  beyond  the 
design  discharge  by  the  cylindrical  piers  (tests  4  and  6). 


a.  Looking  across  ICS  towards  model  trees  on  floodplain.  h.  Looking  upstream  across  floodplain. 


Figure  8.  Typical  ice  jam  formed  in  model. 
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Figure  9.  Measured  downstream  overturning  moment  (test  5, 
Ml),  showing  peak  impulsive  and  sustained  moments. 


m  uv 

a.  Versus  ice  thickness. 


c.  Versus  gap  between  piers. 


b.  Versus  ice  flexural  strength. 


Figure  10.  Scatter  plots  of  peak  overturning  moments.  The 
measured  moments  are  essentially  independent  of  these  parame¬ 
ters  within  the  ranges  tested. 

Figure  9  shows  a  time  series  record  of  the  downstream 
overturning  moment  measured  on  a  pier  during  ice-jam  for¬ 
mation.  Peak  impulsive  moments  (duration  less  than  1  sec¬ 
ond)  were  usually  caused  by  impacts  of  individual  ice  floes, 
while  peak  sustained  loads  developed  during  unsteady  ice- 
jam  formation  and  collapse  events  (which  occurred  on  the 
order  of  10-100  seconds).  Montgomery  et  al.  (1980)  and 
Montgomery  and  Lipsett  (1980)  provide  convincing  evidence 
that  even  massive  bridge  piers  can  respond  to  impulsive  loads 
caused  by  ice  impacts.  The  natural  frequency  of  a  pier  (typ¬ 
ically  1-20  Hz)  governs  its  response  to  these  dynamic  loads 
(Montgomery  et  al.  1980).  Because  the  ICS  piers  will  carry 
no  superstructure,  they  should  have  frequencies  at  least  as 
high  as  bridge  piers,  and  even  ice-impact  loads  will  effec- 
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Figure  11.  Cumulative  probability  for  measured  peak  down¬ 
stream  overturning  moments,  The  straight  line  represents 
the  best-fit  normal  distribution  to  the  data. 


lively  be  long  in  duration.  Thus,  we  selected 
the  largest  downstream  moment,  impact  or  sus¬ 
tained,  measured  on  any  of  the  five  instrument¬ 
ed  piers  as  the  peak  downstream  moment  for  each 
test,  Md. 

Table  2  shows  the  resulting  values  of  M^. 
Interestingly,  these  peak  moments  are  not  cor¬ 
related  with  ice  thickness,  ice  strength,  or  pier 
gap  (Fig.  10).  We  also  saw  no  dependence  on 
pier  location.  This  allows  us  to  treat  the 
values  as  a  single  population  and  thereby 
develop  a  statistical  basis  for  design  loads.  Fig¬ 
ure  11  shows  the  resulting  probability  distri¬ 
bution  for  Md  plotted  on  a  normal-probability 
scale.  A  straight  line  fits  the  data  reasonably 
well,  implying  a  normal  distribution. 

We  also  compiled  for  each  test  the  peak 
downstream  force,  Fx,  peak  transverse  over¬ 
turning  moment,  Mx,  and  the  downstream  over¬ 


t 

Figure  13.  Effective  moment  arm,  Lp,  versus  downstream 
moment,  My,  at  the  time  of  the  peak  downstream  force,  Fx, 
during  each  test. 


turning  moments  at  the  times  of  those  peaks, 

My,  all  measured  by  the  six- axis  load  cell.  As  with  Md, 
Fx  and  Mx  are  not  correlated  with  ice  thickness,  ice 


strength,  or  pier  gap. 
From  the  ratio  My/Fx, 
we  may  compute  the 
effective  moment  arm, 
Zp,  relating  the  mea¬ 
sured  downstream 
force  to  the  down¬ 
stream  overturning 
moment  (Fig.  12). 
This  calculation  as¬ 
sumes  that  vertical  ice 
forces  contribute  neg¬ 
ligibly  to  My,  a  reason¬ 


My 


Riverbed 


Figure  12.  Effective  moment 
arm,  Lp,  relating  downstream 
ice  force,  Fx,  to  downstream 
overturning  moment,  My, 
measured  on  six-axis  load  cell. 


able  assumption  for  cylindrical  piers  during  breakup. 
Figure  13  shows  that  Lp  tends  to  increase  slightly  for 
increasing  My,  indicating  that  rising  ice-contact  height 
and  increasing  downstream  force  both  contribute  to 
large  downstream  overturning  moments.  Similarly,  Fig¬ 
ure  14  shows  that  the  ratio  M^/My  tends  to  decrease  for 
increasing  My,  although  the  correlation  is  not  strong. 

In  only  one  case  (test  3)  did  the  ICS  not  arrest  the 
initial  breakup  ice  run.  This  test  began  with  no  ice  in 
the  main  channel  upstream  of  the  ICS.  The  test  out¬ 
come  agrees  with  observations  at  the  Hardwick  sloped- 
block  ICS:  after  release  of  an  initial  jam,  a  subsequent 
ice  run  consisting  of  1-ft-thick  floes  (or  less)  will  not 
jam  across  the  14-ft  gaps  at  the  ICS.  We  repeated  this 
test  with  thicker  ice  (tests  9  and  10)  and  foxmd  that  the 
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Downstream  Moment,  M  <M  f  t-(b) 
y 

Figure  14.  Ratio  of  transverse-to-downstream  moments  versus 
downstream  moment,  Wy,  at  the  time  of  peak  transverse  moment, 
Mx,  during  each  test. 


cylindrical-pier  ICS  arrested  the  initial  runs  but  released 
the  resulting  ice  jams  at  about  50%  of  the  design  dis¬ 
charge.  However,  this  scenario  should  not  occur  at  the 
prototype  ICS.  The  piers  will  form  a  small  pool  that 
will  preferentially  freeze  early  and  solidly.  With  this 
ice  upstream,  the  ICS  will  arrest  an  ice  run  and  hold  a 
jam  well  beyond  the  design  discharge.  Any  later  ice 
runs  will  impact  and  add  to  the  jam  rather  than  inter¬ 
acting  directly  with  the  piers. 

Design  ice  loads 

The  model  data  can  provide  design  ice  loads  for 
downstream  and  transverse  overturning  moments  and 
forces.  The  peak  downstream  moments  measured  during 
each  test  were  cast  in  a  statistical  framework  for  the 
purpose  of  determining  these  design  loads.  We  recom¬ 
mend  designing  for  loads  at  a  cumulative-probability 
level  of  99.9%  and  describe  design  loads  here  based  on 
that  probability.  This  corresponds  to  an  average  return 
interval  of  1000  years,  assuming  one  ice-breakup  event 
per  year.  Design  loads  based  on  higher  or  lower  proba¬ 
bilities  can  be  determined  using  the  same  procedure. 

The  results  on  Figure  11  directlyyieldMd=2.0x  10^ 
ft-lb  as  the  design  value  for  downstream  overturning 
moment.  We  may  use  this  value  to  estimate  the  other 
design  loads.  Because  Lp  tends  to  increase  slightly  for 
increasing  My,  we  may  obtain  a  conservative  estimate 
of  design  value  of  the  downstream  force,  Fx  =  M^/Lp, 
by  using  the  average  value  of  -  4.4  ft.  This  yields 
Fx  =  0.45  X  10^  lb  as  the  design  value  for  downstream 
force.  Similarly,  because  M^fMy  tends  to  decrease 
for  increasing  My,  we  may  obtain  a  conservative  esti¬ 
mate  of  design  value  of  the  transverse  moment, 
Mx = (Mx/My)  •  Md,  by  using  the  average  value  of  M^IMy 
-  0.45.  This  yields  Mx  =  0.90  x  10^  ft-lb  as  the  design 
value  for  transverse  moment.  Roughly  speaking,  the 


design  downstream  force  acts  halfway  up  the  pier,  and 
the  design  transverse  moment  is  half  of  the  design  down¬ 
stream  moment. 

These  design  loads  may  occur  simultaneously  on  a 
pier  because  they  often  derive  from  the  same  event. 
For  this  reason,  although  we  did  not  record  transverse 
forces,  Fy,  we  may  estimate  the  design  value  as 
Fy  =  Mx/Ip  =  0.20  X  10^  lb.  This  assumes  that  trans¬ 
verse  and  downstream  moments  act  at  the  same  height 
on  the  pier  (4.4  ft  above  the  bed)  for  the  design  event. 

We  may  compare  these  results  with  the  downstream 
loads  used  to  design  the  original  weir-with-piers  ICS. 
Recall  that  its  function  was  different.  During  breakup, 
floes  arriving  from  upstream  would  collect  in  the  pool 
upstream  of  the  ICS,  causing  the  overlying  ice  sheet  to 
ride  up  the  piers,  essentially  to  their  tops.  The  design 
value  of  the  downstream  ice  load  was  assumed  to  result 
from  the  crushing  failure  of  a  1-ft-thick  sheet  across 
the  3-ft  width  of  the  pier  at  a  pressure  of 262  psi,  result¬ 
ing  in  Fx  =  0. 1 1  X  10^  lb  (U.S.  Army  1986b).  Although 
this  is  lower  than  Fx  determined  here,  it  acted  much 
higher  on  the  pier  so  that  the  design  overturning  moment 
was  similar,  M^  =  1.5  x  10^  ft-lb.  Rather  than  experi¬ 
encing  ice-crushing  failure,  however,  the  piers  in  the 
present  ICS  are  subjected  to  direct  ice-floe  impacts  and 
unsteady  ice-jam  formation  and  collapse  events. 

It  is  also  interesting  to  compare  the  recommended 
design  loads  with  those  resulting  from  application  of 
the  AASHTO  (1998)  design  standards  for  bridge  piers. 
These  standards  rely  heavily  on  the  work  of  Montgom¬ 
ery  et  al.  (1 984),  which  also  formed  the  basis  for  design 
standards  in  Canada  (CSA  1988). 

The  AASHTO  standards  differentiate  among  crush¬ 
ing  (across  the  full  width  of  a  pier),  bending,  and  im¬ 
pact  ice-failure  modes.  They  note  that  impact  failure, 
where  a  small  floe  comes  to  rest  before  it  has  crushed 
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over  the  full  width  of  the  pier,  can  be  the  controlling 
mode  on  small  streams  (less  than  300  ft  wide).  This 
mode  provides  the  closest  analogy  to  the  ice-arrest  func¬ 
tion  of  the  proposed  ICS.  For  vertical  piers,  the  appli¬ 
cable  equation  is  (AASHTO  1998) 

Fi  =  C^’P't'W’Ki  (1) 

where 

Ca  =  (5f/w+ 1)0-5  (2) 

and  p  =  effective  crushing  strength 
t  =  ice  thickness 
w  =  pier  width 
K\  =  impact  reduction  factor. 

Here,  w  =  5.0  ft,  and  we  may  assume  t  =  2.0ft  for 
design  purposes.  Equation  2  then  yields  Ca  =  1 .7.  The 
largest  value  of  p  suggested  in  the  standard  is  about 
220  psi,  recommended  for  cases  where  breakup  can 
occur  when  the  average  ice  temperature  is  measurably 
below  the  melting  point  (e.g.,  mid-winter  breakup). 
Lastly,  K{  depends  on  Alr^^  where  A  is  the  floe  plan  area 
and  r  is  the  radius  of  the  pier  nose.  The  model  tests  sug¬ 
gest  that  floes  larger  than  about  30  ft  across  seldom  strike 
a  pier  at  high  velocity.  For  this  case,  K\  ~  0.7  applies. 
These  values  combine,  via  eq  1,  to  predict  F{  =  0.4  x  10^ 
lb,  remarkably  close  to  the  design  value  suggested  here. 

The  AASHTO  standards  do  not  provide  quantita¬ 
tive  guidance  on  the  elevation  at  which  to  apply  the 
ice-impact  force  to  the  pier,  especially  as  these  stan¬ 
dards  are  not  intended  for  the  design  of  ice-control  struc¬ 
tures.  For  this  purpose,  the  model  tests  probably  pro¬ 
vide  the  best  guidance,  Lp  =  4.4  ft.  Also,  because  they 
apply  to  bridge  piers,  the  standards  recommend  simul¬ 
taneous  application  of  the  downstream  force  and  a  trans¬ 
verse  force  equal  to  only  15%  of  that  value.  The  simul¬ 
taneous  transverse  force  recommended  here,  based  on 
model  data,  is  45%  of  the  downstream  force  and  is  con¬ 
servative  compared  to  the  standards. 

Additional  ICS  features 

The  cylindrical-pier  ICS  uses  the  adjacent  treed 
floodplain,  located  along  the  right  bank  of  Cazenovia 
Creek  (see  Fig.  6),  as  a  flow-bypass  channel.  This  pre¬ 
sents  additional  design  requirements:  preventing  sedi¬ 
ment  scour  along  the  right  bank,  preventing  floodplain 
flow  from  reentering  the  main  channel  upstream  of  the 
ICS,  and  retaining  the  ice  pieces  in  the  main  channel  as 
flow  diverts  onto  the  floodplain.  We  may  satisfy  these 
requirements  with  additional  ICS  features. 

Right-bank  riprap 

Model  tests  of  the  cylindrical-pier  ICS  and  the 
sloped-block  ICS  (Lever  et  al.  1997)  revealed  that  flood¬ 


plain  flow  could  reenter  the  main  channel  just  upstream 
of  the  ICS,  where  water  levels  are  below  bank-full.  This 
flow  can  erode  the  ice  floes  arched  between  the  flood- 
plain  and  the  nearest  structural  element,  causing  early 
release  of  the  ice  jam.  Furthermore,  this  flow  is  suffi¬ 
ciently  fast  to  scour  natural  bank  materials,  threatening 
catastrophic  release  of  the  jam. 

Riprap  placed  along  the  right  bank  near  the  ICS  can 
prevent  both  problems.  The  riprap  should  extend  at  least 
one  river  width  (about  1 50  ft)  upstream  and  downstream 
of  the  ICS  and  be  tied-off  to  existing  ground.  To  con¬ 
trol  reentering  flow,  the  top  of  the  riprap  should  be  about 
2  ft  above  the  local  floodplain  elevation,  making  it  EL 
643  here.  A  similar  riprap  emplacement  at  the  Hard¬ 
wick  ICS  has  worked  well  to  prevent  scour  and  reen¬ 
tering  floodplain  flow. 

Although  dynamic  ice-jam  surges  can  occur  near  the 
ICS,  the  banks  do  not  experience  persistent  ice  action. 
Consequently,  local  flow  velocities  rather  than  ice  action 
should  govern  the  riprap  design.  At  Hardwick,  the  riprap 
consists  of  granite  stone  with  D^q  =  1 .7  ft,  and  it  has 
shown  no  signs  of  failure.  For  the  Cazenovia  Creek  ICS, 
D50  =  2.0  ft  should  be  adequate. 

ice-retaining  posts 

The  treed  floodplain  adjacent  to  the  ICS  allows  flow 
to  bypass  the  structure  but  retains  ice  in  the  main  chan¬ 
nel.  However,  several  gaps  exists  in  the  trees  along  the 
right  bank  at  the  Cazenovia  Creek  ICS  site.  In  these 
gaps,  wooden  posts  should  be  installed  to  retain  ice 
pieces  in  the  main  channel. 

Posts  installed  in  the  physical  model  (see  Fig.  6) 
retained  most  ice  floes  in  the  main  channel  during  sim¬ 
ulated  breakup  events.  Maximum  ice  elevations  along 
this  right-bank  “tree  line”  were  5-7  ft  above  the  local 
floodplain  elevation.  We  did  not,  however,  measure 
loads  on  the  posts.  Much  of  the  ice  was  grounded  along 
the  bank,  and  the  water  level  drop  across  the  posts  was 
much  less  than  5  ft. 

Wooden  posts  placed  6  ft  on-center  along  the  gaps 
in  the  right-bank  trees  are  adequate  to  retain  ice.  The 
posts  should  protrude  4  ft  above  the  top  of  the  bank. 
We  may  conservatively  assume  that  a  5-ft  hydrostatic 
head  acts  across  each  post  and  that  the  resulting  force 
acts  2  ft  above  the  ground.  Thus,  each  post  must  resist 
12,000  ft-lb  of  overturning  moment. 

Dry,  select-grade  Douglas  Fir  and  Southern  Pine 
(long  leaf)  have  allowable  tensile  stresses  for  static  loads 
of  about  1800  psi  (CRC  1970).  Thus,  10-in. -diameter 
posts  of  these  woods  are  adequate  to  resist  the  estimat¬ 
ed  overturning  moment;  use  of  12-in.-diameter  posts  is 
conservative.  Note  that  allowable  stresses  may  be  in¬ 
creased  by  100%  for  impact  loads  (CRC  1970),  which 
should  accommodate  moments  from  ice  impacts. 
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ICS  design  recommendations 

We  recommend  an  ICS  consisting  of  nine  5-ft- 
diameter  x  10-ft-tall  cylindrical  piers  spaced  with  12-ft 
gaps  to  reduce  ice-jam  flood  damages  along  Cazenovia 
Creek.  The  ICS  should  be  located  at  the  same  cross 
section  as  the  original  weir-with-pier  ICS,  and  the  top 
of  the  cylindrical  piers  should  be  at  EL  644  ft. 

The  design  ice  loads  on  each  pier  (at  a  cumulative- 
probability  level  of  99.9%)  can  act  simultaneously  and 
are  as  follows: 

= 2.0  X 10^  ft-lb  (downstream  overturning  moment) 
Fx  =  0.45  X  10^  lb  (downstream  force) 

Mx  =  0.90  X  10^  ft-lb  (transverse  overturning  moment) 
Fy  =  0.20  X  10^  lb  (transverse  force). 

Riprap,  consisting  of  stone  with  of  about  2.0  ft, 
should  be  installed  along  the  right  bank  at  least  150  ft 
upstream  and  downstream  of  the  ICS  to  prevent  scour 
and  ice  erosion.  The  top  of  the  riprap  should  be  at  EL 
643  ft. 

Dry,  select-grade  12-in. -diameter  wooden  posts 
should  be  installed  6-ft  on-center  along  the  gaps  in  the 
trees  along  the  right  bank  near  the  ICS.  These  posts 
should  protrude  4  ft  above  the  local  top-of-bank  eleva¬ 
tions.  Insofar  as  possible,  existing  trees  on  the  flood- 
plain  should  not  be  disturbed  during  construction. 

UPSTREAM  EFFECTS 
Upstream  issues 

The  cylindrical-pier  ICS  will  arrest  ice  arriving  from 
several  miles  upstream.  It  will  retain  the  resulting  ice 
jam  at  discharges  much  higher  than  the  maximum  dis¬ 
charge  at  which  an  ice  jam  would  remain  in  the  natural 
channel.  While  this  should  significantly  reduce  ice-jam 
flood  damages  downstream,  water  levels  during  ice 
events  will  be  higher  upstream  of  the  ICS  than  under 


existing  conditions.  The  project  could  thus  affect 
upstream  property  owners. 

Cazenovia  Creek  flows  through  a  steep-walled  val¬ 
ley  from  East  Aurora  to  Mill  Road.  There  are  few  struc¬ 
tures  located  within  this  valley  that  are  vulnerable  to 
flooding.  Table  3  lists  such  structures  from  the  ICS  to 
Transit  Road,  the  upstream  limit  of  this  analysis.  Fig¬ 
ure  15  shows  the  topography  of  this  reach.  Most  prop¬ 
erties  that  abut  the  creek  do  not  have  structures  within 
the  valley. 

The  objectives  here  are  to  determine  the  maximum 
water  level  caused  by  the  ICS  for  all  locations  upstream 
and  to  determine  the  location  beyond  which  the  ICS 
essentially  has  “no  effect.”  We  have  taken  this  latter 
condition  to  be  the  location  where  the  maximum  water 
level  induced  by  the  ICS  is  below  the  100-year  open- 
water  elevation  under  existing  conditions.  The  infor¬ 
mation  sought  here  is  needed  to  determine  the  real  estate 
requirements  (e.g.,  flow  easements,  buy-outs,  flood 
protection)  to  construct  the  ICS.  We  used  a  numerical 
ice-hydraulic  model,  with  input  from  the  physical 
model,  to  determine  upstream  effects  of  the  ICS. 

Ice  breakup  with  ice-control  structure 

The  physical  model  tests  indicated  that  the  ICS  will 
arrest  a  breakup  ice  run  and  retain  the  resulting  ice  jam 
to  discharges  exceeding  6000  cfs  without  catastrophic 
ice  releases  (see  Table  2).  For  the  recommended  ICS 
(12-ft-gaps),  1-ft-thick  ice  of  moderate  strength  should 
start  to  wash  out  through  the  gaps  at  7000-8000  cfs. 
The  ICS  should  retain  stronger  or  thicker  ice  to  higher 
discharges.  For  example,  the  12-ft-gap  ICS  held  1.0-ft- 
thick  strong  ice  to  8600  cfs,  even  after  we  manually 
released  ice  at  one  of  the  gaps  (test  21).  Also,  the  10-ft- 
gap  ICS  held  1. 1-ft-thick  strong  ice  and  1.5-ft-thick 
moderate-strength  ice  to  10,000  and  1 1,000  cfs,  respec¬ 
tively,  without  major  releases  (tests  16  and  19).  Smaller 
floes  did  wash  out  slowly  (through  the  ICS  gaps  and 
onto  the  floodplain)  during  all  tests,  reducing  the  ice 


Table  3.  Structures  within  the  Cazenovia  Creek  valley  from  the 
ICS  site  to  Transit  Road.  Most  properties  that  abut  the  creek  do 
not  have  structures  within  the  valley. 

Distance  Flood 
upstream  of  elevation 

Name _ Descnption _ ICS  (miles)  (ft-NGVD) 


Kotecki  Grove 

Commercia!,  picnic  areas 

0.27-0.36 

648-650 

Leydecker  House 

Single  residence 

1.1 

660 

Leydecker  Road 

Bridge 

1.2 

666 

(low  steel) 

Winspear  Subdivision 

Several  residences 

1. 8-2.0 

672 

Brady’s  House 

Single  residence 

2.1 

676 

Upstream  House 

Single  residence 

2.2 

676 
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Figure  15.  Topographic  map  of  reach  of  Cazenovia  Creek  valley  modeled  for  upstream  effects  caused 
by  the  cylindrical-pier  ICS.  (From  USGS  Orchard  Park  Quadrangle,  1965.) 


volume  upstream  of  the  ICS.  Nevertheless,  we  would 
expect  the  ICS  to  retain  thick,  strong  ice  well  beyond 
the  6000-cfs  design  discharge. 

The  ICS  jarh  will  be  grounded  near  the  structure  and 
free-floating  beyond  about  200  ft  upstream.  Beyond  the 
grounded  section,  the  ICS  jam  is  hydraulically  similar 
to  a  natural  jam,  adding  a  thick,  rough,  floating  ice  cover 
to  the  river.  The  jam  produces  much  higher  water  lev¬ 
els  along  its  length  than  would  occur  for  the  same  dis¬ 
charge  during  open-water  events.  However,  washouts 
and  melting  of  smaller  floes  will  reduce  the  retained 
ice  volume  (and  hence  reduce  jam  length)  throughout 
the  event.  Beyond  about  8000-10,000  cfs,  the  washout 
rate  should  accelerate,  depending  on  the  ice  thickness 
and  strength. 

Ice  jam  volume 

The  volume  of  ice  in  the  jam  retained  by  the  ICS 
governs  the  extent  of  high  water  induced  upstream.  This 
ice  volume  will  change  during  a  breakup  event,  and  we 
require  an  estimate  of  the  volume  retained  as  a  func¬ 
tion  of  discharge,  )^(0.  The  initial  ice  run  from  up¬ 
stream  will  leave  ice  floes  along  the  banks  (in  the  form 


of  shear  walls)  and  on  the  small  floodplains  that  line 
the  valley  floor.  Once  the  ICS  forms  a  jam,  melting  and 
slow  washouts  will  reduce  its  volume  throughout  the 
event.  At  high  discharge,  ice  washouts  will  accelerate 
as  floes  break  and  release  through  the  ICS  gaps.  We 
require  estimates  of  each  term  in  this  process  to  deter¬ 
mine  Vj(Q). 

River  ice  supply 

Buffalo  District  personnel  conducted  ice-thickness 
measurements  in  1983  and  aerial  surveys  in  1983  and 
1985  to  observe  the  ice-formation  process  and  deter¬ 
mine  the  ice  supply  in  Cazenovia  Creek  (U.S.  Army 
1986b).  They  estimated  that  the  pre-breakup  ice  sup¬ 
ply  “from  the  headwaters  to  Cazenovia  Street”  was  7.8 
X  10^  in  1983  and  9.2  x  10^  ft3  in  1985.  Although 
the  ice-thickness  data  themselves  are  not  available, 
maps  used  for  the  surveys  during  1983  record  about 
0.7-ft-thick  solid  ice.  Consequently,  the  thickness  of 
solid  ice  in  1985  was  probably  about  0.9  ft.  These  vol¬ 
ume  estimates  are  consistent  with  ice  covering  100% 
of  the  16  miles  of  river  from  Cazenovia  Street  to  East 
Aurora,  assuming  an  average  cover  width  of  120  ft. 


The  river  length  from  the  ICS  site  to  East  Aurora  is 
about  11  miles.  Along  this  reach,  the  average  cover 
width  should  be  about  100  ft.  A  good  estimate  of  the 
maximum  average  ice  thickness  immediately  prior  to 
breakup  is  1.5  ft.  This  yields  an  ice  volume  of  9  x  10^ 
ft^,  which  we  may  increase  to  10  x  10^  ft^  to  allow  for 
minor  ice  supply  from  the  east  and  west  branches  (U.S. 
Army  1986b).  This  becomes  our  best  estimate  for  the 
pre-breakup  ice  supply,  Vj-. 


Transport  losses 

Field  measurements  have  shown  that  the  ice  volume 
in  a  breakup  jam  can  be  substantially  less  than  the  vol¬ 
ume  of  the  pre-breakup  ice  cover  on  the  contributing 
river  reach.  The  formation  of  shear  walls  and  the  strand¬ 
ing  of  ice  floes  on  small  floodplains  are  the  main  trans¬ 
port  losses  occurring  during  an  ice  run.  The  transport 
loss  coefficient,  Q,  is  defined  as  the  ratio  of  ice  vol¬ 
ume  lost  during  a  breakup  run  to  the  volume  of  the  pre¬ 
breakup  ice  cover.  It  is  thought  to  increase  with  the 
length  of  contributing  reach  (U.S.  Army  1982). 

Table  4  lists  values  of  the  Q  available  in  the  litera¬ 
ture  and  those  calculated  here  for  the  1972  and  1985 
ice  jams  on  Cazenovia  Creek.  We  derived  the  latter 
values  using  reported  jam  locations  and  lengths  (U.S. 
Army  1972,  Predmore  1985).  The  transport  loss  coef¬ 
ficient  is  related  to  the  ice-supply  and  ice-jam  parame¬ 
ters  as  follows: 


B;  h;  .  , 


(3) 


where  L  -  length  of  the  jam 
B  =  width  of  the  jam 
hj  =  thickness  of  the  jam 
p  =  porosity  of  the  jam,  assumed  to  be  0.4 
(Prowse  1990) 

Lr  =  length  of  the  contributing  ice  cover 
-  width  of  the  contributing  ice  cover 
t{  =  average  ice  thickness  of  the  contributing 
ice  cover. 


For  ice  jams  on  smaller  rivers,  typical  values  of  Zj/Lf  ~ 
0. 1 ,  B'JBx  ^  2,  h'ylti  ~  4  result  in  Q  ~  0.5. 

Nearly  all  of  the  available  estimates  of  Q,  including 
ones  based  on  jams  in  Cazenovia  Creek,  equal  or  exceed 
0.3.  The  exception  is  for  an  ice  jam  on  the  Winooski 
River  (which  flooded  Montpelier,  Vermont,  in  March 
1992),  whose  contributing  reach  was  very  short  and  had 
few  overbank  areas  to  strand  ice.  Therefore,  although 
an  average  value  of  0.5  could  be  justified,  we  will  use 
Q  =  0.3  as  a  conservative  estimate  for  transport  losses. 

Ice  melting 

Ice  breakup  normally  occurs  as  a  result  of  rapid 
snowmelt  or  rainfall  during  periods  of  above-freezing 
air  temperatures.  Upstream  of  an  ice  jam,  air-water 
convection,  solar  radiation,  groundwater  inflow,  fluid 
friction,  and  geothermal  flux  will  all  add  heat  to  the 
river.  Under  the  assumption  that  little  ice  remains 
upstream  or  that  it  has  been  stranded  above  the  water 
level  during  the  breakup  surge,  this  heat  flux  will  raise 
the  temperature  of  water  entering  the  head  of  the  jam. 
The  resulting  melting  (and  weakening)  of  ice  floes  within 
a  jam  can  be  significant  during  long-duration,  high- 
discharge  events.  Although  loss  of  jam  volume  may  not 
be  visibly  obvious,  these  effects  probably  couple  with 
rising  hydrodynamic  forces  to  release  natural  ice  jams. 

Table  5  presents  measurements  of  water  tempera¬ 
ture  entering  breakup  ice  jams.  Because  of  the  high  flow 
velocities  and  ice-jam  roughness,  nearly  all  the  avail¬ 
able  heat  is  transferred  to  ice  melting  within  about 
1  mile.  If  the  jam  is  longer  than  this,  ice  melting  at  the 
upstream  end  will  predominate.  As  the  jam  shortens, 
some  warm  water  will  persist  to  the  toe  and  melt  or 
weaken  ice  preferentially  along  main  flow  paths.  In  the 
case  of  a  jam  at  an  ICS,  this  process  contributes  to  its 
ultimate  release  or  washout. 

The  sensible  heat  of  water  is  1  Btu/lbm  ®F  or  about 
62  Btu/ft^  °F.  The  latent  heat  of  fusion  for  ice  is  144 
Btu/lbm  or  about  8200  Btu/ft^  for  solid  ice  (specific 
gravity  0.92).  Assuming  100%  of  the  sensible  heat  of  the 


Table  4.  Transport  loss  coefficients  for  breakup  ice  jams. 


Contributing 

Transport 

reach, 

loss  coefficient 

Reference 

River 

(mi.) 

C( 

Comments 

Calkins  (1978) 

Ottaquechee  R. 

26 

0.9 

Assumed  jam  porosity  of  0.4 

14 

0.9 

(Prowse  1986). 

Prowse  (1986) 

Liard  R. 

300 

0.8 

Entire  event. 

94 

0.4 

Last  24  hours  of  movement. 

Cumming-Cockburn  (1986) 

Credit  R. 

9 

0,5 

3  years  of  field  surveys. 

Tuthill  et  al.  (1996) 

Winooski  R. 

3 

0 

Short  reach,  few  overbank  areas. 

This  work 

Cazenovia  Cr. 

12 

0.3-0.5 

1972  ice  jam. 

16 

0.3-O.5 

1985  ice  jam. 

14 
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Table  5,  Measurements  of  water  temperature  entering  breakup  ice  jams.  The  heat-transfer  length  is  the 
distance  from  the  head  of  a  jam  to  the  point  where  the  water  has  lost  >  90%  of  its  sensible  heat. 

Entering  water  Heat-transfer 
temp.,  AT  length 


Reference 

River 

rF-32) 

(miles) 

Comments 

Calkins  (1984) 

Ottaquechee  R. 

1.3 

0.8 

Upstream  of  refrozen  jam.  time  between 
breakup  and  measurement  unknown. 

Prowse  and  Marsh  (1989) 

Liard  R. 

3.1 

2 

Measured  during  breakup  event. 

Beltaos  etal.  (1998) 

Matapedia  R. 

4.5 

0.2 

Time  between  breakup  and  measurement 
unknown. 

This  work 

Cazenovia  Cr. 

1.8 

unknown 

Measured  by  USGS  about  12  hours  after 
peak  of  1985  ice-jam  hydrograph. 

water  entering  a  jam  is  lost  to  melting  ice,  we  relate  the 
volumetric  melt  rate  of  ice,  Vm,  to  discharge,  Q,  via: 

0.008  *g(cfs)  .ArCF)  (4) 

where  AT  is  the  water  temperature  difference  above 
32°F.  Thus,  water  entering  a  jam  1°F  above  freezing 
will  cause  a  volumetric  melt  rate  of  about  1%  of  river 
discharge.  The  data  in  Table  4  support  a  temperature 
difference  at  least  this  high.  The  expected  rise  in  AT  as 
the  event  proceeds  should  compensate  for  loss  of  sen¬ 
sible  heat  through  the  toe  of  the  jam  as  it  becomes  short¬ 
er.  Therefore,  we  will  use  Fjj^(cfs)  =  0.01  •  as  our 
best  estimate  for  the  volumetric  melt  rate  of  the  ice  jam 
retained  by  the  ICS. 

Ice  washouts 

During  the  model  tests,  we  did  not  quantify  the  rate 
of  ice  loss  attributable  to  major  washouts  or  ultimate 
releases  through  the  ICS.  However,  observations  dur¬ 
ing  washouts  or  releases  at  high  discharge  suggest 
approximately  1%  ice  concentrations  downstream  of 
the  ICS.  Washouts  of  smaller  ice  floes  through  the  ICS 
and  onto  the  floodplain  also  took  place  throughout  the 
tests  without  release  of  the  larger  floes  arched  at  the 
ICS.  For  simplicity,  however,  we  will  assume  that  the 
washout  rate,  is  zero  below  8000  cfs,  and  increases 
to  F^(cfs)  =  0.01  •  g(cfs)for  discharge  above  8000  cfs. 
By  neglecting  washouts  at  low  discharge,  this  approach 
is  probably  conservative  in  its  effect  on  ice  jam  vol¬ 
ume. 

Hydrograph  rise  time 

To  determine  ice-jam  volume  lost  up  to  a  given  dis¬ 
charge,  we  must  integrate  the  loss  rates  (expressed  in 
terms  of  river  discharge)  with  respect  to  time.  This 
requires  an  expression  for  the  rise  time  of  a  character¬ 
istic  hydrograph:  the  slower  the  rate  of  rise,  the  more 
ice  volume  is  lost  to  melting  and  washouts  up  to  a  par¬ 
ticular  discharge. 


We  used  the  hydrographs  from  the  1972  and  1985 
ice-jam  events  (see  Fig.  2)  to  estimate  the  rise  time  for 
discharge  after  ice  jam  formation  (i.e.,  from  just  after 
the  initial  spike  to  the  broad,  ice-free  peak).  This  yielded 
400  cfs/hr  and  500  cfs/hr  for  the  1972  and  1985  events, 
respectively.  We  will  use  a  straight-line  hydrograph  that 
rises  at  500  cfs/hr  from  a  jam-formation  discharge  of 
2000  cfs.  Although  a  hydrograph  rises  more  slowly 
during  the  early  portion  of  an  event,  use  of  the  faster, 
near-peak  rate  produces  a  conservative  estimate  of  the 
ice  volume  lost  to  melting  and  washouts. 

Ice-jam  volume  versus  discharge 

We  may  combine  the  preceding  terms  to  estimate 
the  volume  of  ice  in  the  jam  retained  by  the  ICS  as  a 
function  of  discharge,  lj(0.  Transport  losses  of  30% 
reduce  the  estimated  10  x  10^  ft^  pre-breakup  ice  sup¬ 
ply  to  an  initial  jam  volume  of  Vj  (2000  cfs)  =  7  x  10^ 
ft3.  Following  a  straight-line  hydrograph  that  rises  at 
500  cfs/hr,  melting  losses  occur  at  a  rate  of  1%  of  dis¬ 
charge  throughout  the  event.  Above  8000  cfs,  additional 
losses  ascribable  to  washouts  at  the  ICS  occur  at  a  rate 
of  1%  of  discharge. 

Figure  16  shows  the  resulting  ice-jam  volume  at  the 
ICS  as  a  function  of  river  discharge.  Because  losses 
increase  with  increasing  discharge,  ice  jam  volume 
decreases.  Ice  losses  become  particularly  significant 
above  about  8000  cfs,  as  we  would  expect  from  the 
model  tests,  and  that  just  above  11,000  cfs  we  would 
expect  essentially  all  the  ice  to  have  melted  or  washed 
out  through  the  ICS.  This  estimate  of  Vj(Q)  should  be 
conservative,  and  we  used  it  to  constrain  the  ice-jam 
length  to  predict  upstream  water  levels. 

Numerical  ice-hydraulic  model 

Model  formulation 

We  used  HEC-RAS,  the  Corps’  numerical  hydrau¬ 
lic  model  (U.S.  Army  1998a),  to  calculate  water  sur¬ 
face  profiles  through  the  ICS  reach  for  both  open-water 
and  ice-jam  cases.  Briefly,  HEC-RAS  treats  water  flow 
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Figure  1 6.  Volume  of  ice  in  the  jam  retained  by  the  cyiindrlcal>pier 
iCS  as  a  function  of  river  discharge.  Transport  losses  reduce  10,000,000 
ft^  pre-breakup  ice  supply  to  an  initial  jam  volume  of  7,000,000  ft^  at  2000 
cfs.  Ice  melting,  plus  washouts  through  the  ICS  above  8000  cfs,  reduce 
ice-jam  volume  as  discharge  increases  throughout  an  event. 


as  one-dimensional,  steady,  and  gradually  varied  (spa¬ 
tially),  as  did  its  predecessor  HEC-2.  At  each  cross  sec¬ 
tion,  it  can  include  an  ice  cover  of  known  thickness  or 
solve  for  the  thickness  of  an  ice  jam  according  to  stan¬ 
dard  theory  (U.S.  Army  1998b).  The  latter  is  a  steady- 
state  theory  that  treats  the  ice  as  a  granular  material  with 
no  cohesion.  The  program  solves  the  one-dimensional 
force  balance  for  the  jam,  where  the  under-ice  water 
shear  and  jam  self-weight  are  resisted  by  the  shear 
strength  of  ice  at  the  banks.  The  ice  jam  floats  on  the 
water  and  adds  a  rough  top  surface,  resulting  in  water 
levels  that  are  much  higher  than  for  open  water  at  the 
same  discharge.  The  user  can  constrain  the  jam  to  the 
main  channel  (appropriate  for  treed  floodplains)  and 
can  set  hydraulic  roughness  (Manning’s  n)  of  the  ice  to 
a  fixed  value  or  allow  it  to  vary  with  jam  thickness. 
The  program  calculates  the  volume  of  ice  in  a  jam, 
allowing  us  to  simulate  a  jam  of  known  volume  at  each 
discharge  (upstream  of  which  is  open  water). 

Input  to  the  numerical  model  included  50  cross  sec¬ 
tions  along  3.6  miles  of  Cazenovia  Creek,  from  600  ft 
downstream  of  Mill  Road  to  100  ft  downstream  of  Tran¬ 
sit  Road  (see  Fig.  15).  The  cross  sections  came  from 
three  sources:  1)  an  HEC-2  deck  called  Cazpfe  that  used 
surveyed  cross  sections  from  1984,  2)  closely  spaced 
surveyed  cross  sections  obtained  in  1998  to  construct 
the  physical  model,  and  3)  supplemental  cross  sections 
determined  from  topographic  maps  (1  in.  =  200  ft) 
obtained  from  the  Town  of  West  Seneca.  As  far  as  pos¬ 
sible,  we  checked  these  data  for  consistency  in  terms 
of  vertical  and  horizontal  alignment.  The  HEC-2  deck 
included  the  geometry  of  the  two  bridges  in  the  reach, 


at  Mill  Road  and  Leydecker  Road.  We  input  the  5-ft- 
diameter  xlO-ft-tall  x  12-ft-gap  cylindrical-pier  ICS 
using  the  HEC-RAS  feature  “multiple  blocked  obstruc¬ 
tions.” 

All  model  runs  used  subcritical  flow.  The  down¬ 
stream  rating  curve  derived  from  the  output  of  the  HEC- 
2  deck  (RS  308.00,  U.S.  Army  1986b).  Several  cali¬ 
bration  flows  were  used  to  compare  computed  water- 
surface  elevations  (WSE’s)  with  physical-model  and 
field  data.  Final  runs  were  made  using  the  1 00-year  ( 1  %- 
exceedence)  open-water  flow  of  1 5,600  cfs  and  ice-jam 
flows  of 4000-13,000  cfs. 

Model  calibration 

Open-water  WSE  data  without  the  ICS  were  avail¬ 
able  from  three  sources:  1)  field  measurements  along 
the  ICS  site  in  1984  at  discharges  of  783  and  3170  cfs, 
2)  a  HEC-2  computed  water-surface  profile  at  2000  cfs 
(U.S.  Army  1986b),  and  3)  a  water-surface  profile  at 
1 5,600  cfs  from  a  recent  flood  insurance  study  (FEMA 
1992).  As  shown  in  Figure  17,  minor  adjustment  of 
channel  and  over-bank  roughness  values  along  the  reach 
gave  good  agreement  between  the  calibration  data  and 
the  HEC-RAS  model  results  at  these  four  discharges. 

Open-water  WSE  data  with  the  ICS  were  available 
from  the  physical  model  at  discharges  of  1550,  2900, 
and  5800  cfs.  We  used  these  to  adjust  the  contraction 
and  expansion  coefficients  in  the  HEC-RAS  model  to 
calibrate  the  water-level  change  across  the  ICS  (Fig. 
18). 

No  WSE  field  data  were  available  for  ice-jam  con¬ 
ditions  upstream  of  Mill  Road,  probably  because  this 
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has  not  been  an  area  of  concern  during  past  events.  The 
ice  run  through  the  reach  in  1 985  caused  5-  to  6-ft-high 
shear  walls  upstream  of  both  Mill  Road  and  Leydecker 
Road  (Predmore  1985),  although  these  probably  reflect 
very  dynamic  jam  formation  and  collapse  processes  at 
the  onset  of  the  event  (say,  2000-4000  cfs).  Also,  local 
residents  (J,  Durand  and  J.  Brady)  noted  that  ice  runs 
have  caused  below-flood-stage  elevations  at  and  up¬ 
stream  of  Winspear  Road,  but  they  differed  in  opinion 
about  whether  stationary  jams  formed. 

The  physical  model  provided  WSE  data  with  the  ICS 
for  many  different  ice-jam  conditions.  We  selected  three 
data  sets  (2000,  6000,  and  6100  cfs)  that  offered  fairly 
long,  steady-state  ice  jams  upstream  of  the  ICS.  As  with 


all  the  tests,  however,  the  jams  were  grounded  at  the 
ICS,  a  feature  that  cannot  be  modeled  directly  with 
HEC-RAS.  Instead,  we  matched  the  measured  WSE’s 
by  setting  6-fit-thick  ice  covers  along  the  first  50  ft 
upstream  of  the  ICS  (Fig.  19).  We  selected  a  conserva¬ 
tive  value  of  0.08  for  Manning’s  n  of  the  underside  of 
the  ice  jam  and  used  this  for  all  the  final  ice-jam  runs. 
Runs  conducted  with  variable  Manning’s  n  (i.e.,  n 
dependent  on  jam  thickness)  yielded  similar  results. 
Also,  runs  with  the  ice  jam  toe  at  Mill  Road  converged 
with  these  results  about  1000  ft  upstream  of  the  ICS, 
indicating  that  river  geometry  governs  ice-jam  thick¬ 
ness  except  near  the  structure.  At  ice-jam  discharges  of 
9000-12,000  cfs,  we  thickened  the  jam  toe  to  8  ft  and 
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extended  it  to  300  ft  upstream  of  the  ICS  to  achieve 
stable  calculations  near  the  ICS. 

Results 

Figure  20  shows  the  open-water  profiles  for  the  100- 
year  discharge  (15,600  cfs)  with  and  without  the  ICS. 
The  effect  of  the  ICS  disappears  beyond  about  2500  ft 
upstream.  Although  the  ICS  raises  the  100-year  open- 
water  flood  elevation  at  Kotecki  Grove  by  about  0.6  ft, 
ice  jam  conditions  actually  dictate  maximum  WSE  at 
this  location. 

Figure  21  shows  ice-jam  profiles  at  4000, 7000,  and 
10,000  cfs,  with  ice  volume  limited  according  to  dis¬ 


charge  (see  Fig.  16).  HEC-RAS  requires  a  fixed  ice- 
jam  thickness  at  the  head  of  a  jam,  and  we  used  1,5  ft 
(i.e.,  the  parent  ice  coverthickness)  for  all  calculations. 
Open-water  conditions  prevail  upstream  of  a  jam.  Fig¬ 
ure  21  shows  that  jam  length  decreases  dramatically 
with  increasing  discharge,  a  combination  of  decreas¬ 
ing  ice  volume  and  increasing  jam  thickness. 

We  determined  the  maximum  water  levels  expected 
upstream  of  the  ICS  by  assembling  the  open-water  and 
ice-jam  profiles  for  all  discharges  into  a  single  data  file. 
Figure  22  shows  typical  results.  Where  the  profiles 
cross,  the  discharge  that  establishes  the  maximum  WSE 
changes.  Table  6  provides  the  maximum  WSE,  the  dis- 
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Elevastion  (ft)  Elevastlon  (ft)  Elevastion  (ft) 


Main  Channel  Distance  (ft) 
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Main  Channel  Distance  (ft) 


Table  6.  Maximum  waterlevels  expected  upstream  of  the  ICS,  Note  that  the  ICS  has  no 
effect  beyond  about  9600  ft  upstream,  where  the  100-year  open-water  profile  dictates  maximum 
waterlevels.  The  ICS  is  located  at  the  site  originally  proposed  for  the  weir-with-pier  ICS  (USACE 
1986a,  1986b).  identified  as  RS  337.27  in  Cazpfe  HEC-2  deck. 

Distance  Maximum  water  Governing 
upstream  of  ICS  surface  elevation  discharge  Open  water  or 

(ft) _ (ft-NGVD) _ (c^ _ ice  conditions _ Landmark _ 


15962 

684.74 

15600 

Open-Water 

13062 

674.93 

15600 

Open-Water 

11618 

670.10 

15600 

Open-Water 

Upstream  House 

11178 

667.62 

15600 

Open-Water 

Brady’s  House 

10818 

666.15 

15600 

Open-Water 

10403 

665.79 

15600 

Open-Water 

WInspear  Subdivision 

9692 

663.72 

15600 

Open-Water 

Winspear  Subdivision 

7792 

662.98 

7000 

Ice 

6292 

662.71 

7000 

Ice 

.6192 

662.31 

7000 

Ice 

Leydecker  Rd 

6164 

662.23 

7000 

Ice 

Leydecker  Rd 

5950 

661.86 

7000 

Ice 

Leydecker  House 

5650 

661.23 

7000 

ice 

4037 

656.64 

9000 

Ice 

3727 

656.13 

9000 

Ice 

3252 

655.86 

9000 

Ice 

2782 

655.52 

9000 

Ice 

2247 

655.00 

9000 

Ice 

1892 

654.58 

10000 

Ice 

Kotecki  Grove 

1417 

653.52 

10000 

ice 

KoteckI  Grove 

1080 

651.66 

10000 

Ice 

960 

651.15 

10000 

Ice 

897 

650.98 

10000 

Ice 

840 

650.69 

10000 

Ice 

720 

650.02 

10000 

Ice 

600 

649.73 

11000 

Ice 

480 

649.50 

11000 

Ice 

360 

649.29 

11000 

Ice 

299 

648.96 

11000 

Ice 

240 

648.84 

11000 

Ice 

132 

648.10 

11000 

Ice 

120 

647.98 

11000 

Ice 

50 

647.13 

11000 

Ice 

0 

645.34 

15600 

Open-Water 

ICS 
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charge  that  establishes  that  WSE,  and  governing  sur¬ 
face  condition  (open-water  or  ice-jam)  at  each  cross 
section. 

At  the  Winspear  Subdivision  and  above,  the  100-year 
open-water  profile  dictates  maximum  water  levels.  That 
is,  the  ICS  has  no  effect  on  water  levels  beyond  about 
9600  ft  upstream.  From  there  to  below  the  Leydecker 
House,  the  7000-cfs  ice-jam  profile  dictates  maximum 
water  levels.  The  governing  ice-jam  discharge  then 
increases  with  proximity  to  the  ICS,  such  that  10,000 
cfs  dictates  the  maximum  water  level  expected  at 
Kotecki  Grove.  Of  the  structures  within  the  Cazenovia 
Creek  valley  (Table  3),  only  Kotecki  Grove  and  the 
Leydecker  House  are  affected  by  ice  jams  at  the  ICS. 
Note  also  that  ice  jams  at  the  ICS  do  not  contact  the  bridge 
at  Leydecker  Road  (low-steel  elevation  666.00  ft). 

DISCUSSION 

The  cylindrical-pier  ICS  is  a  refinement  of  the 
sloped-block  ICS  installed  in  Hardwick,  Vermont.  These 
wide-gap  structures  arrest  breakup  ice  runs,  retain  the 
partially  grounded  ice  jams,  and  use  existing  treed  flood- 
plains  to  bypass  flow.  They  function  well  without  requir¬ 
ing  a  weir  or  prepared  floodway.  Cylindrical  piers  with 
12-ft  gaps  should  retain  ice,  particularly  thin  ice,  at  much 
higher  discharge  than  sloped  blocks  with  14-ft  gaps. 
For  breakup  events  that  pose  the  greatest  flood  threat 
(thick,  strong  ice,  heavy  rain,  or  rapid  snowmelt),  the 
cylindrical-pier  ICS  should  retain  ice  until  well  after 
ice  has  cleared  from  downstream  reaches.  Thus,  this 
new  ICS  should  substantially  reduce  ice-jam  flood  dam¬ 
ages  along  Cazenovia  Creek  in  West  Seneca. 

Construction  plans,  specifications,  and  cost  estimates 
were  completed  for  the  original  weir-with-piers  ICS  (US 
Army  1986a).  The  most  expensive  items  were  com¬ 
mon  excavation,  compacted  fill,  and  concrete  with 
forming  (for  the  upstream  pool,  prepared  floodway,  and 
weir-with-piers  structure,  respectively).  The  new  ICS 
eliminates  the  first  two  items  and  substantially  reduces 
the  third.  We  may  estimate  approximately  50%  cost 
savings  for  the  remaining  items  (e.g.,  clearing  and  grub¬ 
bing,  roadways,  riprap,  engineering,  supervision,  land). 
The  construction  cost  for  the  new  ICS  might  thus  be 
about  one-quarter  to  one-third  of  the  cost  of  the  origi¬ 
nal  ICS.  The  new  ICS  will  not  have  a  gate  to  operate 
and  will  not  require  dredging  to  remove  deposited  sedi¬ 
ment.  Despite  these  lower  costs,  model  results  suggest 
that  it  should  perform  at  least  as  well  as  the  original 
concept. 

The  calibrated  HEC-RAS  model  should  provide 
reliable  estimates  of  the  water-surface  profiles  expected 
for  both  open-water  and  ice-jam  events,  with  and  with¬ 
out  the  ICS.  Because  ice-jam  water  levels  are  so  much 


higher  than  even  100-year  open-water  ones,  the  jam 
length  or  volume  essentially  governs  the  upstream  influ¬ 
ence  of  the  structure.  Insofar  as  possible,  we  have  based 
the  terms  that  affect  jam  volume  on  data  from  the  physi¬ 
cal  model,  the  literature,  or  ice-jam  events  in  Cazenovia 
Creek  itself,  and  have  selected  values  at  the  conservative 
end  of  the  ranges  for  each  term.  Consequently,  the  pre¬ 
dicted  maximum  water  levels  caused  by  the  ICS  should 
be  conservative. 

CONCLUSIONS 

We  recommend,  as  the  new  Cazenovia  Creek  ICS,  a 
structure  consisting  of  nine  5-ft-diameter  x  10-ft-tall 
cylindrical  piers  spaced  with  12-ft-gaps  across  the  main 
chaimel.  The  treed  floodplain  should  be  left  intact  to 
act  as  a  flow-bypass  channel.  About  300  ft  of  riprap  is 
needed  along  the  right  bank  near  the  structure  to  pre¬ 
vent  scour  and  to  delay  reentry  of  floodplain  flow  until 
farther  downstream.  Comparison  of  the  model  results 
indicates  that  the  new  ICS  should  have  an  ice-retention 
capability  at  least  as  high  as  the  original  weir-with-piers 
ICS.  For  ice  thickness  approximately  1  ft  or  greater, 
we  would  expect  the  ICS  to  arrest  a  breakup  ice  run 
and  retain  the  ice  jam  for  a  discharge  exceeding  7000 
cfs.  Slow  washouts,  rather  than  catastrophic  releases, 
will  likely  be  the  release  mode  at  higher  discharges  or 
with  thinner  ice.  We  saw  no  performance  advantage 
for  a  cylindrical-pier  ICS  with  10-ft  gaps,  whereas  one 
with  14-ft  gaps  allowed  ice  washouts  more  easily. 
Because  the  new  ICS  does  not  include  a  weir  and  exca¬ 
vated  pool,  and  makes  use  of  the  existing  floodplain  as 
a  bypass  channel,  it  should  be  substantially  cheaper  than 
the  original  weir-with-piers  ICS. 

Interestingly,  the  analysis  for  upstream  effects  has 
revealed  the  ice-retaining  capacity  of  the  ICS  as  a  bal¬ 
ance  between  two  needs:  the  need  to  protect  downstream 
areas  from  natural  ice-jam  flooding  with  the  need  to 
minimize  upstream  flooding  by  jams  at  the  ICS.  Melt¬ 
ing  and  slow  washout  of  ice  through  the  ICS  reduces 
extent  of  flooding  upstream  without  endangering  flood¬ 
ing  downstream.  Within  our  present  capabilities 
for  physical  and  numerical  modeling,  the  12-ft-gap, 
cylindrical-pier  ICS  appears  to  strike  a  good  balance 
of  these  competing  needs. 
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